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X-RAY PROPERTIES OF THE Be/X-RAY SYSTEM 2S0114+650 = LSI+ 6 BOOIO 


a. Koenigsberger^, J.H. Swank, A.E. Szymkowiak^, and N.E. White^ 
Laboratory for High Energy Astrophysics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 20771 

ABSTRACT 

Results are presented from experiments on the Einstein Observatory, 

HEAO-1 and OSO-8 on the temporal and spectral properties of 2S0114+650. In a 
12 hour Einstein MPC and SSS observation two episodes of flaring occurred by 
an order of magnitude over about 1 hour* Variability on shorter t1m.s scales 
showed a preferred period of 14.9 minutes, but periodic pulsations were not 
seen in the HEAO-1 and OSO-8 data. The spectra above 1 keV were typical of 
accreting X-ray pulsars, a power law E"“ with a ~ 1.2 for the number spectrum 
up to a cut-off at ^ 14 keV. There is some evidence that the spectrum is 
steeper when the source is quiescent at low luminosity, but no large spectral 
changes attended the flares. Absorption column densities were- consistent with 
interstellar reddening of the proposed companion, the Be star LSI+65°010. In 
the minute and hour variability and in the spectral character, 2S0114+650 is 
similar to other Be Si^ar-neutron star binary X-ray sources. Variations over 
several days in the OSO-8 data suggest orbital effects. 

^Dept. of Astronomy, Penn State University 
2 a1so Dept. Physics & Astronomy, Uni'tJ. of Maryland 


le INTRODUCTION 

There a class of low luminosity (1033-1035 erg s~M 2-10 keV X-ray 
sources associated with main sequence Be stars (Bradt et al. 1977; Maraschl, 
Treves, and van den Heuvel 1976; Rappaport and van den Heuvel 1981). Stable 
X-ray pulsations have been detected from many of these systems and long term 
measurements of these pulsation periods have Indicated that the X-ray pulsars 
must be In orbits about the Be stars with periods of order tens to hundreds of 
days (McClIntock et al. 1977; White et al. 1978; LI et al. 1979; White et al. 
1982). Evidence for such a period has been obtained only for the transient 
X-ray source 4U0115+63 (Rappaport et al. 1978). Periodic outbursts from some 
of these systems also suggest similar orbital periods (e.g. Watson, Warwick 
and Ricketts 1981). Secular changes in the periods of the pulsars, mass loss 
considerations, and the peak luminosities seen from the transients Indicate 
that the pulsars are accreting neutron stars (White, Mason and Sanford 1976; 
Rappaport et al. 1978; Rappaport and van den Heuvel 1981; White et al. 

1982). Correlated evolution of the pulse periods and the companion mass loss 
has been speculated (van den Heuvel 1977; Ghosh and Lamb 1978). The range of 
binary conditions Imply a range of accretion rates that may Influence the 
beaming and spectra of the pulsars (Rappaport and Joss 1977; White, Swank and 
Holt 1983). Rappaport and van den Heuvel (1981) have suggested on 
evolutionary grounds that neutron star binaries with Be star companions are 
likely to have longer periods than those with OB supergiant companions. Thus 
the properties of sources of this type have a bearing on several Interesting 
questions. 

The X-ray source 2S01 14+650 was first reported in the SAS-3 galactic 
survey by Dower and Kelly (1977). Its optical counterpart. Identified by 
Margon (1977) as LSI+65°010 (Hardorp et al . 1959), is an 11th magnitude 
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B0.5III star which presents a broad Ha emission line In Its spectrum (Margon 
1980). There are two other hard X-ray sources associated with dwarf Be stars 
within 3° of 2S01 14+650; y Cas (MX0053+64) and the recurrent transient 
4U01 15+63 (Bradt, Oaxsey and Jenigan 1979). The distance estimates for y Cas 
of 300 pc (Moffat et al. 1973); for ISOl 15+634 of 3.5 kpc (Rappaport et al . 
1978); and for 2S0114+650 of 1.4 kpc (Margon 1980) put y Cas In the local 
0>^1on arm and the other two In opposite sides of the Perseus spiral arm (e.g. 
Humphreys 1979). Their proximity In the sky reflects a high probability of 
finding Be stars in regions dense In early type stars. A similar clustering 
of X-ray binaries with early companions Is found in Centaurus where we view 
along the Carina spiral arm (White et al. 1980). While the X-ray properties 
of y Cas and 4U01 15+63 are now reasonably well studied (e.g. White et al. 

1982; Rose et al. 1980) those of 2S0114+650 are not. In this paper we report 
the spectral and temporal properties of this source and find them to be very 
similar to other members of this class. 

II. OBSERVATIONS AND RESULTS 

2S01 14+650 was observed by the Einstein Observatory on 1979 Jan 17 for a 
duration of 12 hours with the Solid State Spectrometer (SSS; 0.5-4. 5 keV) and 
the Monitor Proportional Counter (MPC; 1.5-10 keV). Descriptions of these 
Instruments were given by Joyce et al. (1978) and Giacconi et al. (1979). The 
m.oan observed flux in the 1.5-10 keV energy band was 7 MPC cts s”l, 
corresponding to ~ 2 x lO"^^ erg cm"2s~l and, assuming a distance of 1.4 kpc, 
a luminosity of ~ 4 x 1034 ergs s"l. The HEAO-1 A2 experiment (see Rothschild 
et al. 1979) scans during 11-17 August 1978 detected 2S0114+650 at a flux 
level down by a factor of 5 with respect to the MPC observation. Twenty-two 
days earlier, on July 20, 1978, a 6 hour pointed observation of the 3.6 s 
transient X-ray pulsar 4U01 15+63 was made (cf. Rose et al. 1979). Although 
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4U0115+63 was the target, no X-ray pulsations at 3.6 s were detected, and 
2S0114+650 probably provided the counts we detected. The 2-60 keV "A" 
detector of the Goddard Cosmic X-ray Spectroscopy Experiment on OSO-8 
(Serlemitsos et al. 1976) scanned the source during 27 Jan.-l Feb. 1976. 

Again the response during the scans indicated that 4U01 15+63 was not detected 
during the observation, 
a. Temporal Behavior 

Figure 1 is a plot of the MPC 2.56 s rates data binned so that each value 
corresponds to 40.96 seconds. Two prominent outbursts are evident. In each 
case there is an increase in flux by a factor of 3 above the mean and lasting 
for ~ 1 hr. This represents a range in 2-10 keV X-ray luminosities of 4 
X 10^3 - 6 X lO^*^ ergs s"*^. There is also variability on timescales of 
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shows considerable noise at low frequencies that reflects the flaring 
activity. On top of this there is evidence for significant power at 14.9 
min No other Significant pulsations were found for shorter periods, with 
j upi,>er limits of 9% down to 11. seconds, 17% for 11. - 0.04s and 22% for 0.04 - 

1 0.004 s (Weisskopf and Leaky, private communication.) 

In Figure 3 have plotted four portions of the 3-10 keV MPC data folded 
on the 14.9 min period using in each case the same initial epoch. The four 
I . plots are a) first outburst region, b) second outburst, c) prior to first 

outburst, d) between outbursts. {These intervals are indicated in 

* 

Figure 1). Although the amplitudes, mean intensities and shapes of the pulses 


vary (as would be expected from the raw data), all four curves peak between 
phases 0.3 and 0.5 and have a minimum around phase 0.7, indicating rough 
coherence in the pulses over the whole data set. The expected positions of 
pulse maxima in the raw data, based on the profiles in Figure 3a, are marked 
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In Figure 1. Most of the pulses are visible, although their varying structure 
is apparent, especially in those cases in which there are multiple peaks 
within a period. Figure 4(a) shows the pulse profile of 2S01 14+650 obtained 
by folding all the MPO data with the 14.9 minute period, while Figure 4(b) is 
the hardness ratio defined as flux (3^10 keV)/flux(l .5-3 keV). As in X Per 
(White et aU 1982) there is a hardening of the spectrum during decline that 
ends at pulse minimum. Note however, the hardness ratio in 2S01 14+650 is ~ 3 
times larger than that in X Per corresponding to a flatter spectrum, a ~ 1.2 
versus ~ 2 for X-Per in a power law model ,, as we discuss below. The 
peak-to-mean amplitudes for the energy bands 0.5-4. 5, 1. 5-3,0 keV and 3-10 keV 
are 20%, 26% and 21%, respectively. 

The OSO-8 and HEAO-1 scanning observations, each covering 6 days, are 
shown in Figure 5. Variability on timescales of several days can be seen, 
especia'^ly in the OSO-8 light curve where a transition from maximum to minimum 
occurs over 4 days and represents a factor of ~ 10 decrease in flux. There 
is evidence for more rapid variability which is probably similar to the 
flaring activity observed in the Einstein observations (Figure 1). During 
the 1 day interval Jan 31 , 1975 (= JD 2,443,808.5 in Figure 4) the source 
was undetectable with a 3a upper limit of ~ 2.1 x 10”^^ ergs crtT^ g-l. 

Figure 6 shows the light curve for the 6*^ HEAO-1 point, assuming 
4U0115+63 was not contributing. The average flux implied is the same as that 
during the scanning observations and during the low portions of the Einstein 
MPC data. There is no evidence in the data for pulsations at ~ 3.6 s, the 
pulse period of 4U0115+63, although the upper Vimit of 30% is not very 
restrictive. Variations of a factor of 2 with time scales of 10-30 min are 
evident. But there is also no evidence for pulsations at 14.9 minutes greater 
than ~ 20%, compared to the level of ~ 25% seen in the MPC data. We were also 
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unable to confirm the pulsations in the OSO-8 data for which the flux was 
high, although our upper limit of 25% Is comparable to the observed MPC 
modulation, 
b. Spectrum 

Spectral data Indicate that for all the observations the 2-10 keV 
spectrum was hard, similar to that from other members of this class (Becker et 
al. 1979; White et al . 1982; White, Swank and Holt 1983). Table 1 shows the 
fits for the absorbed power law and thermal bremsstrahlung models that fit 
both the MPC 2-10 keV data and the SSS 1.3-4. 5 keV data. There was no 
significant difference between the spectra of the flares and of the low states 
during the Einstein observation. However, the spectrum during the HEAO-1 A2 
point, when the source was weaker, appears vo be slightly steeper. Figure 7 
shows Inferred spectra from these observations for the data above 1.3 keV. 

The OSO-8 data (from the early part of the observation, when the source was 
bright) are only shown above 10 keV. Selow this there was some contribution 
to the OSO-8 pulse height data from the Tycho supernova remnant, which has a 
(2-10 keV) spectrum softer than the spectrum of 2S01 14+650 seen during the 
Einstein observation (Pravdo et al. 1979). We deduce from the observed 
spectrum that the high state power law for 2S01 14+650 has a cutoff at least 
above 15 keV, 

The column densities obtained in the proportional counter fits are ~ 3 . 
1022 cm~2 and are consistent, within 3a, with the reddening of LSI+65°010. 
Margon (1977) reports B-V ~ 1.18, which would Imply column densities ~ 10^^ 
cm"2 (Gorensteln 1975). Because of source confusion in the Tow energy OSO-8 
pulse height data, we cannot determine whether there Is any evidence for 
absorption increases that might attend a possible eclipse transition prior to 
Jan. 31, 1975, when the source was below the OSO-8 detection limit. However, 
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there 1s no strong Indication of enhanced absorption during the HEAO-1 A2 
observation, when again the source was weak. 

There Is some discrepancy In Figure 7 between the lowest energy (1-3 keV) 
SSS and MPC estimates of the Incident spectra. This Is due to the strong 
model -dependence of the lowest energy MPC points (channel 1). Also, +*'ie SSS 
spectra are affected by an excess below 1 keV, Because of uncertainties In 
the background subtraction needed for this observation, it is difficult to 
determine the properties of this soft excess. However, If 1t Is assumed to be 
emission from the primary itself then it corresponds to an ~ 10^3 ergs s”^ 
with a temperature of ~ 10® K. This Is more than three orders of magnitude 
larger than that expected from a dwarf B star (Pall avid n1 et al. 1981). 

III. DISCUSSION 

The X-ray speetfUiii and time variability of 2S0114+650 are very similar to 
other members of the X Per class of Be stars. Two outbursts, separated by ~ 6 
hours, and representing an order of magnitude variation In luminosity, were 
resolved during the 12-hour Einstein Observation. In its flaring activity, 
2S01 14+650 Is similar to X Per where outbursts lasting about 1 hour have been 
known to occur at ~ 11 hour intervals (White et al. 1976), From the OSO-8 
observation of 2S01 14+650 we see that the frequency and/or amplitude of 
flaring varies on longer time scales (1-2 days), which may be related to 
orbital phase. Also, the flux can remain below 1 yJy for as long as a day. 
While this mlant suggest an eclipse, simlliarly low flux levels were also seen 
for several hours during the Einstein observation and the HEAO-1 point. Even 
at this low level variations occur which suggest the compact object Is still 
In view as the X-ray source. Also, the spectrum is not significantly more 
absorbed at any time as would be expected for a grazing eclipse by a star with 
an extended atmosphere. Thus It seems unlikely that the low observed 1n the 
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OSO -8 data Is due to an ecHpse, 

Unitts on possible binary periods for 2S011 4+650 can be made by assuming 
that its X-ray luminosity is a result of spherical aiccretion from the stellar 
wind of the primary as described by Davidson and Ostriker (1973). The 
gravitational energy released by the accreted material, the accretion radius, 
the velocity of the wind, and Kepler's third law give the following relation 
for the orbital period 


P(day) ~ 140 


,3/4 j^3/4 


1-3/4 

**33 


R„ -3/4 9/4 . 

''3 


T' 


- 1/2 


where ? is an efficiency factor for the conversion of gravitational potential 

• 

energy into the observed radiation, is the primary's mass-loss rate in 
10"^ Mq yy'~^ (assumed spherically symmetric)^ L 33 is the X-ray luminosity in 
units of 10^3 ergs s"l, and are the radius and mass of the compact 
object, respectively, Vg is the wind velocity in units of 10 ^ km s"l and Mj is 
the total mass of the binary system. A crude estimate of Vg can be infered 
from the width of the Ha emission line if it is assumed that the emission is 
produced in the radially expanding wind of the Be star. From the first 
profile in Margon's (1980) Figure 3 we estimate AX < 42 A which 
implies Vg < 0.9. (This could be an overestimate if the Ha-emission region 
is much larger than the radius of the orbit). Typical values for the wind 
velocities in X Per and y Cas lie between 500 and 1000 km/s with mass-loss 
rates between 10"^ and 'i0“H Mq yr"^ (Hammerschlag-Hensberge et al . 1980; Snow 
1981). 

Using My = 15 Mq, Mx = 1.4 Mq, V 3 ~ 0.9, L 33 = 10, Rx = 10 km for M = 
lO’ll - 10“8 Mq yr’^, we obtain possible orbital periods in the range ~ (0.07 
- 18) day. If ? ~0.1, only orbital periods < 3 d could give sufficient 
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X-ray luminosity. In contrast to the long periods (>15 d) Indicated for other 
Be star systems. Longer periods would be consistent If V 3 < 0 . 9 , M > 10“® Mq 
y"^, or the mass loss Is concentrated In the orbital plane. A similar problem 
was encountered by White et al. (198?) their calculation of the mass-loss 
rates required to power the observed X-ray luminosities In X Per and y Gas. 

The variations In luminosity In 2S011 4+650 on timescales of days may either 
reflect changes In the mass loss rate from the star, an eccentric orbit for 
the neutron star, or a misalignment between the orbital plane and the 
equatorial plane of the rotating Be star. 

The Einstein data exhibit a preferred period of 14.9 minutes for the 
short term variability with the pulsations maintaining coherence over the 44 
cycles contained 1n the observation. This period Is close to that of the 13.9 
min period of the prototype X Per. We were not able to confirm periodicity In 
the other data although the limits on pulsations were not very restrictive. 
This periodicity is reminiscent of another Be star X-ray source with a neutron 
star companion, y Cas, White et al. (1982) report variations which look like 
pulses and appear coherent over short Intervals but over multiple observing 
Intervals change period. They suggest that these variations result from 
Inhomogenitles 1n the wind. 

The spectrum of 2S0114+650 1s similar, as measured by i.he MPC and SSS, to 
typical pulsar spectra above 1 keV (White, Swank and Holt 1982). The cutoff 
at > 15.6 keV indicated by OSO-8 when the spectrum is flat 1s also typical. 

For example, 4U1 145-61 and 41)1538-52 exhibit cutoffs ^;t 16 keV. Although 

the power law index derived from the Einstein and OSO-8 data Is somewhat 
smaller than those of X-Per and y Cas, the HEAO A2 spectrum would be typical 
of the latter. The long term change between the Einstein and HEAO A2 
observation Is the first marked spectral change reported for these systems. 
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The Tack of spectral correlation with relatively shot't time-scale flares (15 
min) exhibited in the Einstein data is also typical of the behaviors of X-Per 
and y Cas (White et al. 1982). While on the basis of y Cas, X-Per, 4U1 145-61, 
and 4U1 538-52 it would be tempting to conclude that the steepness of the power 
law reflects the much lower luminosities of y Cas and X-Per, 2S01 14+650 
provides a counter-example. Indeed tho observations suggest the power law for 
low luminosity states is not unique. 

IV. CONCLUSIONS 

We have shown that 2S01 14+650 shares several characteristics with other 
members of the Be/X-ray class. In particular, its temporal behavior is 
analogous to that of X-Per and y Cas, although its spectrum is somewhat 
flatter. The variability on 14.9 minutes timescale in the Einstein data 
suggests the possibility that 2S01 14+650 contains a pulsar, as in the case of 
X-Per. However, further observations are needed to determine whether this 
period persists over longer timescales. 

If the binary period is > 15 days, as is suggested by other Be/X-ray 
sources, then mass-loss rates > 10"® yr"^ or material concentrated in the 
orbital plane are required to produce the observed X-ray luminosity. The 
X-ray luminosity shows a wide spectrum of variability time scales from the 15 
minute pulsations to several days. It is not known whether the long term 
changes are periodic and thus helpful in determining the orbital parameters, 
or whether they reflect changes or inhomogeneity in the mass-loss of the Be 
star. 
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FIGURE CAPTIONS 

Figure 1 - Raw MPC (1.5-10 KeV) light curve from day 382.87 to 383.37 of 
1978 (JD 2443891). The intensity scale is given in MPC counts/2.56 sec 
(1 MPC ct s“^ « 1.6 X 10"11 epg cm~2 s“^. The arrows point to the 
expected positions of maxima according to the 14.9 minute period and are 
labeled with the corresponding cycle number. a,b,c,d are intervals 
referred to in Figure 3. 

Figure 2 - Fourier transform of the MPC data for periods longer than 2 
minutes. 

Figure 3 - Pulse profiles obtained by folding the 3-10 keV MPC data with P = 
893.8 seconds and initial (arbitrary) epoch at Day 382.87, 1978. a) 
first outburst, b) second outburst, c) prior to the first outburst and d) 
between outbursts. 

Figure 4 - a) Pulse profile from the entire, folded data set in Figure 1. 
b) Hardness ratio. 

Figure 5 - Light curves during the OSO-8 observation and the HEAO-1 A2 

scan. For the OSO-8 intensities 1 count/s is approximately 1 yjy. This 
corresponds approximately to 0.0025 cts cm"^ s“^* for the HEAO-1 A2 rate 
combination used. The scales have been chosen so that they are 
approximately equivalent. The OSO-8 points are averages over ~ 2000 s, 
the A2 points over 0.1 day. 



16 


Figure 6 - Light curve for the HEAO-1 A2 pointed observation. 

Figure 7 - Inferred Incident spectra above 1,3 keV for all observations. 

The times covered In the MFC data Include the Intervals of SSS data, but 
are not Identical. That the OSO-8 spectrum matches on to the MFC 
spectrum shown Is presumably fortultlous, since both are averages over 
apparently random flares. 
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TABLE 1; SPECTRAL FITS 


Thermal Brems^ Power Law^ Fj^{2-10 keV) L^^ 


Detector 

kT(keV) NH(1022cm-2) 

a Nnno22cm“2) (I0"10ergs 

cm"2s‘"l) (103^ ergs s"M 

MPC 

Flare 

>25 

4.5 

1.2±0.1 

3.2*1 3.6 

22.9 

(1.3-12 keV) 

Low 

>20 

4.5 

1.2±0.3 

3.2tl 0.2 

1.3 

SSS 

Flare 

>10 

3±2 

1.0±0.5 

2±1 


(1.3-4. 5 keV) 




1.2(Fixed) 

1.6±0.2 



Low 

>10 

3±1 

0.7±0.1 

0.3±0.2 






1 .2(Fixed) 

0.5±0.2 


HEAO-1 A2 

Low 

14±9 

3±2 

1 .8±0.2 

5±2 



(2-20 keV) 

OSO-8 High 40±13 — l.2(F1xed)C 

(10-60 keV) 


®(dn/dE)j = c g E"^ where g Is the Gaunt factor. cr(E) is the cross 

section for photoionization of cold gas (Fireman 1974). 


^ {dn/dE)p * cE"“ 

(l for E < E^ 

^ dn/dE = (dn/dE)p x- -(E-E^)/E- 

Le ^ ^ for E > E^ 

Estimato of Jtot^ 1 umi nosi ty assuming the spectrum c for the hard component with 
and Ep from the OSO-8 results. 
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